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Abstract 

Normal-state magnetic properties of cuprate high-T c superconductors are inter- 
preted based on the self-consistent solution of the t-J model of Cu-0 planes. The 
solution method retains the rotation symmetry of spin components in the paramag- 
netic state and has no preset magnetic ordering. The obtained solution is homoge- 
neous. The calculated temperature and concentration dependencies of the magnetic 
susceptibility are close to those observed in experiment. These results offer explana- 
tions for the observed scaling of the static uniform susceptibility and for the changes 
in the spin correlation length, spin-lattice and spin-echo decay rates in terms of the 
temperature and doping variations in the spin excitation spectrum. 
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Magnetic properties of cuprate perovskites have been extensively studied dur- 
ing the last years, both because of their unusual behavior and in the hope that 
they might provide insight into the physical origin of high-T c superconductiv- 
ity [1,2,3,4,5,6,7]. Considerable progress has been made in this field with the 
use of phenomenological approaches, exact diagonalization of small clusters 
and - for heavily overdoped materials - in a RPA treatment. However, the 
basic issues of the magnetic behavior of underdoped cuprates, where strong 
electron correlations reveal themselves in full measure, have not yet been com- 
pletely resolved. In this paper magnetic properties of underdoped cuprates in 
the normal state are investigated using the t-J model of Cu-0 planes and the 
method of Ref. [8] for calculating Green's functions. This method has merits of 
retaining the rotation symmetry of spin components in the paramagnetic state 
and of the absence of any predefined magnetic ordering. Test calculations with 
this method for small clusters and for the undoped case demonstrated good 
agreement with the exact diagonalization and Monte Carlo results. For the 
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considered 20x20 lattice and parameters of cuprates the calculated magnetic 
properties of the model appear to be close to those observed in these crys- 
tals. The results give an insight into mechanisms responsible for the unusual 
magnetic behavior of cuprate perovskites. 

The Hamiltonian of the 2D t-J model reads [9] 

H = E *nmaL<W + \ E J nm + ^m) , W 

nmcr ^ nm 



where a ncr = \na)(n0\ is the hole annihilation operator, n and m label sites 
of the square lattice, a = ±1 is the spin projection, |ner) and |n0) are site 
states corresponding to the absence and presence of a hole on the site. For 
nearest neighbor interactions t nm = — tJ2 a ^n,m+a an d J nm = JJ2 a ^n,m+a 
where t and J are the hopping and exchange constants and the four vectors a 
connect nearest neighbor sites. The spin-| operators can be written as s z n = 
\ J2a cr\na)(na\ and s n = |ner)(n, —<r\. 

The method suggested in Ref. [8] is based on Mori's projection operator 
technique [10] which allows one to represent Green's functions in the form 
of continued fractions and gives a way for calculating their elements. The 
residual term of the fraction can be approximated by the decoupling which 
reduces this many-particle Green's function to a product of simpler func- 
tions. In this way we obtained the following self-energy equations for the hole 
G{kt) = -i0(t)({a ka (t),al a }} and spin D(kt) = -i9(t)([s^(t), si k ]) Green's 
functions: 



n „ x [4Jq(A + 1 + 7k )]- 1 n(M + 4Jd( 7 k ~ 1) 

D(M " ^_n(M-^ ' 

(2) 

G(ku) = <f>[uj - e k + H - S(kcj)] -1 , 

where 7 k = \ X) a ex P(^ka), /i is the chemical potential, = |(1 + x), x is the 
hole concentration and 



wj[ = 16J 2 a|d|(l - 7 k) (A + 1 + 7k), 

(3) 

e k = (40t + 6Ci0 -1 £ - ZF^- 1 J) 7k . 



The parameter of vertex correction a, which improves the decoupling in the 
residual term, is set equal to its value in the undoped case, a = 1.802. The 
parameter A which describes a gap in the spectrum of spin excitations at 



2 



(ir, 7r) [see Eq. (3)] is determined by the constraint of zero site magnetization 
(sf ) = in the paramagnetic state. The constraint can be written in the form 



du coth 



to 
2T 



B(ku), 



(4) 



where B{kuo) = — 7T~ l lmD(kuj) is the spin spectral function, N is the number 
of sites and T is the temperature. Notice that in the considered 2D system the 
long-range antiferromagnetic ordering is destroyed at any nonzero T [11] and, 
as can be seen from the above formulas, at T = and x > 0.02. The value of 
x and the nearest neighbor correlations C\ = (s^s^) and Fi = (a\ai +a ) are 
defined from G(ku) and D(ku) in the usual way. 

The self-energies in Eq. (2) read 



Imn(kw) = ^^(A + 1 + 7 k) E(7k - 7k+k') 2 
iV k' 



x 



J du'[n F {u + J) - n F (u')]A(k + k',cu + J)A(k!J), 



ImS(kcu) = 1 ^ J du)'[n B (-u)') + n F (u - a/)] 



N(p 



(5) 



x 



7k-k' + 7k + sgn(w')(7k-k' - 7k) 1 



l + 7k' 
1 - 7 k / 



1] 1 and A(ku) 



xA(k-k',LU-Lu')B(k'u'), 

where n F (uj) = [exp(u;/T) + l] -1 , ub{oj) = [exp(a;/T) 
— 7r _1 Im G(kuo) is the hole spectral functions. The source of damping of spin 
excitations described by Eq. (5) is the decay into two fermions. Another source 
of damping, multiple spin excitation scattering, is considered phenomenologi- 
cally by adding the small artificial broadening —2r]uj^ i] = 0.02t to Imll(kc<j). 
The broadenings — i] is also added to ImS(kcj) to widen narrow lines and to 
stabilize the iteration procedure. 

The same derivation for the transversal Green's function gives ((sj^lsk" 1 )) = 
2D(ku) indicating that the used approach retains properly the rotation sym- 
metry of spin components in the paramagnetic state. 

For low x and T the bandwidth of the dispersion is small in comparison 
with St, the bandwidth of uncorrelated electrons. This is a manifestation of 
the band narrowing in the antiferromagnetic surrounding. 



Equations (2)-(5) form a closed set which was solved by iteration for the 
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Fig. 1. The dispersion of spin excitations. Vertical bars show decay widths 
|lmli(ku;)|/(2(x>k)- Points Y, M and S correspond to k = (0,7r), (tt,tv) and 
(7r/2,7r/2), respectively. 

parameters t = 0.5 eV, J = 0.1 eV corresponding to cuprates [12]. 
The frequencies of spin excitations satisfy the equation 

uj 2 - RelT(kcj) - ul = (6) 

[see Eq. (2)]. For low x and T their dispersion is close to the dispersion of 
spin waves (see Fig. la). The main difference is the spin gap at the 
magnitude of which grows with x and T (Fig. lb). In an infinite crystal this 
gap is directly connected with the spin correlation length £. Indeed, for large 
distances and low T we find 



<*f *5> = E e * kI / ** coth (^) ^(M 

k 

oce^/llD^e-l 1 !/*, 



(7) 



where Q = (n, n) and ^ = o/(2yA) with the intersite distance a. For low x 
we found that A 0.2x and consequently £ ~ a/y/x. This relation has been 
experimentally observed in La 2 -a;Sr x Cu04 [13]. 

As seen from Fig. lb, with growing x the spin excitation branch is destroyed 
in some region around the V point - Eq. (6) has no solution for real uj due to 
large negative ReII(ku;). For fixed T the size of this region grows with x. In 
the considered model with rise of T the branch is recovered in this region. This 
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Fig. 2. Spin correlations along the diagonal [i.e., 1 = (m,m)] of the crystal for 
x = 0.12. The respective temperatures are indicated near the curves. 

is a consequence of the temperature broadening of the quasiparticle peak in 
the hole spectrum. With this broadening |ReIl| becomes smaller and Eq. (6) 
has again a real solution. 

An example of the variation of the spin correlations C mn = (s%Sq), 1 = (to, n) 
with distance is given in Fig. 2. For large enough x and T the correlations de- 
cay exponentially with distance in the considered finite lattice. As mentioned, 
the method used has no preset magnetic ordering. The character of the or- 
dering is determined in the course of the calculations. As seen from Fig. 2, 
only the short-range antiferromagnetic ordering was found in our calculations. 
Stripes or other types of the phase separation were not revealed. Conceivably 
such phase separations are not connected with the strong electron correlations 
described by the t-J model. 

The magnetic susceptibility is connected with the spin Green's function (2) by 
the relation x z (k^) = — 4:/j, 2 B D(kw), where fj, B is the Bohr magneton. Exper- 
iments on inelastic neutron scattering give information on the susceptibility 
which can be directly compared with the calculated results. Such comparison 
is carried out in Fig. 3. YE^CuaCv-y is a bilayer crystal and the symmetry 
allows one to divide the susceptibility into odd and even parts. For the an- 
tiferromagnetic intrabilayer coupling the odd part can be compared with the 
calculated results. The oxygen deficiencies y = 0.5 and 0.17 in the experimen- 
tal data in Fig. 3 correspond to the hole concentrations x = 0.05 and 0.11, 
respectively [14]. As seen from Fig. 3, the calculated data reproduce correctly 
the frequency dependence of the susceptibility, the values of the frequency 
for which Imx(7r, 7r) reaches maximum and their evolution with doping. The 
growth of the frequency of the maximum with x reflects the respective increase 
of the spin gap. We notice also that the calculated temperature variation of 
the susceptibility is in good agreement with experiment. In absolute units the 
calculated maxima of Inr\;(7r, tt) are 1.5 — 2 times larger than the experimen- 
tal values which is connected with some difference in decay widths of spin 
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Fig. 3. The imaginary part of the spin susceptibility for k = (n, ir). Curves demon- 
strate calculated results for T = 0.02i 116 K, x = 0.043 (a) and 0.08 (b). Squares 
are experimental results obtained in normal-state YBa2Cu307_ y at T = 100 K for 
y = 0.5 (a) and 0.17 (b) [7]. 
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Fig. 4. The uniform static spin susceptibility vs. T. 

excitations. 



The temperature dependence of the uniform static spin susceptibility Xo = 
x(k — > 0,u = 0) is shown in Fig. 4. The calculated values lie in the range 2- 
2.6 eV _1 which is close to the values 1.9-2.6 eV" 1 obtained for YBa 2 Cu 3 7 _. i , 
[5]. The dependence Xo(T) has a maximum and its temperature T m grows with 
decreasing x. Analogous behavior is observed in cuprates for large enough x 
[1,3]. In Fig. 4 T m ps 600 K which is close to T m observed in La2- x Sr x Cu04 for 
comparable x [1]. As known, in the undoped antiferromagnet T m pa J [15]. On 
the high-temperature side Xo(T) tends to the Curie- Weiss dependence l/T. 
The decrease of Xo below T m is sometimes considered as the manifestation of 
the spin gap. In our opinion this statement is incorrect. For moderate x and 
T the long-wavelength part of the spin excitation spectrum does not feel the 
gap at (7r,7r). For small but finite values of k x(k, 0) oc J^dcu' 'B(kw') / 'uj' '. 
The function B(k.cu') has a maximum which is shifted to lower frequencies and 
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loses its intensity with increasing T for such wave vectors. In the above integral 
the maximum is superimposed with the decreasing function 1/u/ which finally 
leads to the nonmonotonic behavior of xo(T)- 

As seen from Fig. 4, the two curves for the different x are very close in shape 
and can be superposed by scaling to the same values of maximum Xo an d 
T m . Analogous scaling was observed in L^-xSr^CuC^ [1]. As follows from 
the above discussion, the source of this scaling is that holes and temperature 
fluctuations lead in a similar manner to the softening of the maximum in 
B{kuo') for long wavelengths. 

The spin-lattice relaxation and spin-echo decay rates were calculated with the 
use of the equations [5] 



1 


1 


a T ip T 


2/4 AT 


1 


0.69 


63^2 
I 2G 


128/4 



k u 

{^E 63 ^(k)[Rex(kO)f (8) 

i-E 63 ^(k)Re X (kO) 
JV k 

where the indices a and (3 in the form factors a Fp{\s) [5] indicate the nucleus 
type and the direction of the applied static magnetic field H, respectively. 
a = 63 corresponds to Cu. The form factor 63 F e is the filter for the Cu spin- 
echo decay time 63 T 2G . Our calculated results and the respective experimental 
data are given in Fig. 5. The calculations reproduce satisfactorily main pe- 
culiarities of the temperature dependencies of the spin-lattice and spin-echo 
decay rates. The growth of ^TiT)~ l with decreasing x is connected with the 
increase of the spectral intensity of spin excitations near (ir, re) which make 
the main contribution to this rate. For the same x ( 63 TiT) _1 is one order of 
magnitude larger than (^T^T) -1 , the spin-lattice relaxation rate at O sites 
(not shown here). This is a consequence of Im% which is strongly peaked near 
and the form factors which test different k regions [5]. The calculated 
spin-lattice relaxation rates are smaller than the experimental values due to 
the approximation made in the calculation of D(ku) which somewhat under- 
estimates Im% at low frequencies. 

For moderate x with increasing T the low-frequency region of Im%(k pa Q), 
Q = (ir, 7r) first grows due to the temperature broadening of the maximum in 
its frequency dependence and then decreases due to the temperature growth 
of the spin gap. It is the reason for the nonmonotonic behavior of ( 63 TiT) _1 in 
Fig. 5a. For moderate x and low T the magnitude of the spin gap is determined 
by the hole concentration and does not depend on T. This temperature range 
corresponds to the growth stage in Fig. 5a. The independence of the gap from 
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Fig. 5. The temperature dependencies of the spin-lattice relaxation and spin-echo 
decay rates at Cu sites. Open circles with right axes represent experimental re- 
sults, filled circles with left axes are our calculations. (a,c) calculations for H||c and 
x = 0.12, measurements in YBa2Cu30g.63[3] (x ~ 0.1 [14]). (b) calculations for 
nonoriented configuration with x = 0.043, measurements in Lai.96Sro.o4Cu04 [2]. 

T means that in this temperature range also £ does not depend on tempera- 
ture which is a distinctive feature of the quantum disordered regime [4,5]. For 
temperatures above the maximum of (^TiT) -1 we found that 63 T 1 T/ 63 T 2 g- ~ 
const (see Fig. 5c) and oc y/A oc T which points to the quantum criti- 
cal z = 1 regime. These results are in agreement with the phenomenological 
treatment of experiment in YB^Cu^Or-y carried out in Ref. [5] and the tem- 
perature of the maximum in Fig. 5a is close to of that work. 

For small x the spin gap grows with temperature starting from low T and 
( 63 TiT) _1 decreases monotonously, as seen in Fig. 5b. Due to the form- factor 
the region near (ir,ir) does not contribute to ( 17 7YT') -1 . There is a cardinal 
difference between the behavior of Im% for k Q and away from (tt, it). Due 
to the spin gap in the former case the frequency of the maximum in Im x(o;) 
increases with T, while in the latter mentioned above, it decreases. 

This frequency softening leads to the growth of the low- frequency Im% and 
(17^7^-1 a t l ow T with their saturation for higher temperatures. Analogous 
behavior is observed in experiment [3]. 

This work was supported by the ESF grant No. 4022. 
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